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■ ABSTRACT 

Q: 

^NJ ' We present a semi-analytic study of the equilibrium models of close binary systems 

' containing a Newtonian fluid star (mass m and radius Ro) and a Kerr black hole (mass 

J> ■ M) in circular orbit. We consider the limit M 3> m where spacetime is described 

by the Kerr metric. The tidally deformed star is approximated by an ellipsoid, and 
£T) . satisfies the polytropic equation of state. The models also include fluid motion in the 

stellar interior, allowing binary models with nonsynchronized stellar spin (as expected 
| for coalescing neutron star-black hole binaries) to be constructed. Our relativistic, 

compressible Roche- Riemann model is a generalization of the incompressible, corotating 
Roche model studied earlier by Fishbone (1973). 

Tidal disruption occurs at orbital radius r ti d e ~ R${M /m) 1 ^ , but the dimensionless 
ratio rtide = ^tide/[-Ro(-^/ m ) 1 ^ 3 ] depends on the spin parameter of the black hole as 
well as on the equation of state and the internal rotation of the star. We find that the 
general relativistic tidal field disrupts the star at a larger rtide than the Newtonian tide; 
. £h \ the difference is particularly prominent if the disruption occurs in the vicinity of the 

black hole's horizon. In general, ftide is smaller for a (prograde rotating) Kerr black 
d ' hole than for a Schwarzschild black hole. 

We apply our results to coalescing black hole-neutron star and black hole-white 
dwarf binaries. The tidal disruption limit is important for characterizing the expected 
gravitational wave signals and is relevant for determining the energetics of gamma ray 
bursts which may result from such disruption. 

Subject headings: stars: neutron - stars: black hole - stars: white dwarf - gravitation 
- hydrodynamics - gamma ray burst 
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1. Introduction 

The problem of tidal disruption of fluid stars by black holes in close binaries was first studied 
in the 1970s (e.g., Fishbone 1973; Mashhoon 1975). Lattimer & Schramm (1976) specifically con- 
sidered neutron star disruption, and focused on mass ejection in such an event and its implications 
for nucleosynthesis. In recent years, coalescing binary systems containing a black hole (BH) and a 
neutron star (NS) or a white dwarf (WD) have become increasingly important both because they 
are one of the most promising sources of gravitational waves (GW) to be detected by ground-based 
laser interferometric GW detectors (LIGO/VIRGO/TAMA/GEO) or spaceborne interferometers 
such as LISA (see Thorne 1997 for a review), and because they potentially serve as a central engine 
for gamma ray bursts (e.g., Narayan, Paczynski & Piran 1992; Meszaros, Rees & Wijers 1998; Piran 
1998). 

Current numerical simulations of NS-BH binaries and WD-BH binaries (Kluzniak & Lee 1998; 
Fryer et al. 1998; Uryu & Eriguchi 1999; Janka et al. 1999) are based on Newtonian theory (in 
some cases, the gravitational field of the BH is approximated by a pseudo-Newtonian potential), 
and mostly (with the exception of Janka et al. 1999) neglect gravitational radiation reaction on the 
dynamics. While such simulations provide useful insight, full understanding of the coalescence will 
have to wait for the development of numerical relativity. 

In this paper, we develop a semi-analytic model of close binaries containing a fluid star (mass m) 
and a massive BH (mass M). To make the problem tractable, we make a series of approximations. 
The obvious simplification is that in the limit M S> m, the spacetime can be approximated by 
the Kerr metric. Furthermore, we assume that the radius of the companion star, Rq, is much 
greater than Gm/c 2 , allowing the self gravity to be treated in Newtonian theory. The companion 
star is modeled using the standard ellipsoid ansatz (Chandrasekhar 1969) and its compressible 
generalization (Lai, Rasio & Shapiro 1993,1994; hereafter LRS1 and LRS2). This has the advantage 
that it reduces the infinite degrees freedom and the partial differential equations to a set of coupled 
ordinary differential equations. These equations are then reduced to coupled algebraic equations by 
demanding quasi-equilibrium. Fishbone (1973) has constructed "relativistic Roche models" where 
the incompressible fluid body is tidally locked with the orbital motion. However, it is necessary to 
consider models where internal fluid motion is allowed, since in a coalescing binary the NS spin is not 
expected to be synchronized with the orbit (Kochanek 1992; Bildsten & Cutler 1992). Newtonian 
Roche- Riemann binary model has been studied in LRS1. The "relativistic Roche- Riemann model" 
(the subject of this paper) represents a generalization to these previous works by allowing for 
nonsynchronized stellar spin and including general relativistic tidal field. We note that Shibata 
(1996) has studied some aspects of the relativistic Roche- Riemann model where the fluid body 
is assumed incompressible. Post-Newtonian correction in binary ellipsoid models has also been 
studied (e.g., Taniguchi & Shibata 1997; Shibata & Taniguchi 1997; Taniguchi 1999). 

Our relativistic Roche-Riemann models (and other semi-analytic models based on ellipsoid) 
clearly have limitations (see also §7). For example, at small orbital radii, higher-order tidal fields 
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are expected to be important and the star will not have an exact ellipsoid shape. In the case of 
NS-BH binaries, the conditions M/m 3> 1 and Rqc 1 / (Gm) 3> 1 are satisfied only approximately. 
However, given that fully relativistic numerical calculations are still prohibitively difficult, our 
analytic models have the virtue in that they provide a wide array of semi-quantitative information 
about the behavior of the these systems. For example, some of the questions we address in this 
paper are: How does the tidal limit in general relativity compare with that in Newtonian gravity? 
How does the tidal limit depend on the spin parameter of the BH (and other parameters related to 
the NS equation of state and internal rotation)? The answers to these questions are important for 
determining the gravitational wave signals and the energetics of the gamma ray bursts that may 
result from the tidal disruption. While the absolute value of the tidal limit calculated from our 
model may not be accurate, its dependence on various parameters should be robust (see §7), and 
can be used to calibrate results from numerical simulations when they become available. 

In §2 and §3 we present the basic equations of our model and derive the scaling relations. 
Numerical results are given in §4, where we focus on the tidal limits for different models. The tidal 
effect on the orbital dynamics is considered in §5, and in §6 we discuss possible applications of 
our results to gamma ray bursts. In §7 we comment on the validity and limitations of our models. 
Throughout the paper we shall use units such that G = c = 1. 



2. Basic Equations 
2.1. Tidal tensor 

The formalism of calculating tidal tensor in general relativity has been studied by many authors 
(e.g., Fishbone 1973; Mashhoon 1975; Marck 1983). The basic calculations and results are briefly 
summarized as follows. Consider a body (with mass m much smaller than the BH mass M) on a 
geodesic in the Kerr metric, (i) One starts with the equation of geodesic deviation: 



d 2 n a 



+ flijAV = o, (i) 



where u is the 4-velocity of the center-of-mass, n is the 4-separation vector, Rp^g is the Riemann 
tensor, which can be calculated in certain convenient basis; (ii) Construct a parallel-transporting 
orthonormal basis carried by the body: {e( ), e^), e( 2 ), e (3)}> with e( ) = u, and satisfies the parallel 
transport equation V u e( a ) = (a = 0, 1, 2, 3); (hi) The geodetic deviation equation can be reduced 
to 

^ + £ C «*i = 0, (1 = 1,2,3), (2) 
j 

where Xi = e^n* 1 , and the tidal tensor is 

Cij = Ra^se^e^eJ^e^y (3) 
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We set up the parallel-transporting frame as in Figure 1, with e( 3 ) perpendicular to the orbital 
plane. The nonzero components of the tidal tensor in this frame is given by (Marck 1983) 



C22 



M 



M 



r 2 + K 9 . 
1 - 3 5— cos 2 , 



1-3 



r 2 + K 



sin 2 * , 



M 



C33 = -^1 + 3^ 



K 



C12 — C21 



M ( 3r 2 + K 
"2 r 2 



sin 2^ 



Here the angle ^ is specified by the equation 

E(L Z -aE)+a 



dr 



r 2 + K 



(4) 
(5) 
(6) 
(7) 

(8) 



where E is the the energy (per unit mass), L z is the z-orbital angular momentum (per unit mass), 
K = (L z — aE) 2 is the Carter constant (with Q = 0), and the dot ("■") denotes d/dr in this and 
all the following equations. In circular orbit, we have (e.g., Misner, Thorne, & Wheeler 1973) 



E 



r 2 - 2Mr + ay/Mr 
r\[P 

VMr(r 2 - 2a\TMr + a 2 ) 



r\[~P 



(9) 
(10) 



where P = r 2 — 3Mr + 2a\jMr. Thus we have 

K = (aE-L z ) 2 , 



(11) 



Note that in general relativity, the angle ^ is not equal to 9, the true anomaly of the orbit, which 
satisfies 



2M 



-aE+\\- 



2M\ 



r J 



L 7 



(r 2 +a 2 -2Mr) . 



(12) 



2.2. Fluid Equations 

Including internal pressure and self-gravity, the Euler equation in the parallel-transporting 
frame reads 

dVi_ ldP d$ v 

where Vi = dXi/dr and <E> is the potential for the (Newtonian) self-gravity. It is convenient to set 
up another frame, tied to the principal axis of the body (see Fig. 1). This "principal frame" is 
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related to the parallel-transporting frame by a rotation of angle <f>: 



(14) 



where Xh = T22 = cos0, T\i = —T21 = sin0, T33 = 1, and the other elements are zero. In the 
principal frame, the Euler equation becomes: 



x + 2fixx + Ox(fixx) + Oxx = --VP - V$ - V$ t idc, 



where 



^tide — r. y ] 



ij 



The tidal tensor in the principal frame is 



The nonzero components are: 



CH = q" 



M 
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C22 = 
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r 2 +K 



cos 2 (^ - 4>) 



sin 2 (^ - 0) 



M 



C33 
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Cl2 = C 2 1 = -^T 



K 



3r 2 + K . 



2 r 2 



sin2(^ - 4>) 



(15) 
(16) 

(17) 



(18) 
(19) 
(20) 
(21) 



2.3. Fluid Equations in the Ellipsoid Representation 

Since the tidal force is a linear function of Xj, for an incompressible fluid, ellipsoidal figures 
describe exactly the solution of the fluid equation (15) (Chandrasekhar 1969). We can further 
generalize the ellipsoid figures to compressible fluid by assuming that the surfaces of constant 
density inside the star form self-similar ellipsoids and that the the velocity of a fluid element is a 
linear function of X{ (LRS1,LRS2). These approximations reduce the infinite degrees of freedom and 
partial differential equations, governing the dynamics of the companion, to five dynamic variables 
associated with the body and a set of ordinary differential equations. The five dynamic variables are 
the three principal axes a±, 02, 03 and two angles 0, A, defined such that dcfi/dr = £1 and dX/dr = A. 
Here f2 is the angular velocity of the ellipsoidal figure measured in the parallel transporting frame, 
and A characterizes the internal fluid motion (see eq. [24] below) via A = — a\a2(,/{ a \ + o, 2 ) with 
C the (uniform) vorticity along the z-axis in the frame corotating with the figure. The physical 
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significance of (j) is illustrated in Fig. 1, but A is without meaning and only its derivative, A, appears 
in our equations. The orbit is specified by the orbital separation r and the true anomaly 6. 

Several parameters must be given to complete the description of the system: In addition to 
the stellar mass m, the BH mass M (satisfying M 3> m) and the Kerr parameter a, there are 
parameters associated with the equation of state and the initial stellar rotation. We adopt the 
polytropic equation of state, 

p = Kp 1+1 / n , (22) 

with n the polytropic index. Note that K (not to be confused with that K in the tidal tensor; see 
eqs. [4]- [7]) is uniquely determined from n, m and R , the radius of the nonrotating star with the 
same mass. Thus instead of K, we shall use Ro as a parameter. Finally, we allow the star to have 
some initial rotation, Q s , when the binary separation is infinite (see §4.3). 

The equations governing the dynamics of the system can be obtained using the variational 
principle. It is convenient to write the Lagrangian as the sum of two terms: 

L = ^body + ^tide- (23) 

We assume that the star moves on a circular geodesic around the BH, and neglect tidal effect on the 
orbital motion — This is an excellent approximation for M S> m even if the star is highly distorted 
(see §5). Thus the orbital terms do not appear in the Lagrangian. The internal Lagrangian of the 
body, Lbodyj is exactly the same as in the Newtonian treatment. The velocity field u of the fluid 
inside the star relative to its center of mass can be written as 



u 



ai A- tt ) x 2 ei + ( -— A + Q ) xie 2 



/ Oi 0,2 (Z3 > 

+ — xiei H x 2 e 2 H x 3 e 3 , (24) 

ai a 2 a 3 



/ V °i 

where ei, e 2 and e 3 are the basis unit vectors along the principal axes of the ellipsoid. On 
integration, the kinetic energy of the star relative to its center of mass is given by 

12 1 
T = -/(A 2 + ft 2 ) - -K n maia 2 Aft + —K n m{a\ + a\ + d 2 ), (25) 

with / = K n m{a\ + a 2 )/5 and where n n is a constant depending on n and specifies the mass 
concentration inside the star (see Table 1 in LRS1). The internal energy is: 

U = Jdmn^- = k x Kp]J n m (26) 

where k\ is another constant depending only on n, and p c oc m/(aia 2 a 3 ) is the central density. The 
self gravitational potential energy is given by: 

3 m? I 

W = 2& ' With J = + ^ 2 ° 2 + A ^ (27) 

where R = (aia 2 a 3 ) 1 ^ 3 is the mean radius of the ellipsoid, and the dimensionless index symbols Ai 
are defined as in Chandrasekhar (1969; §17). The Lagrangian for the star is then 

£bod y = T — U — W. (28) 
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Using the tidal tensor derived earlier (see eq. [16]), the tidal Lagrangian may be calculated as 

Aide = -Wtide = - j d 3 i/)(x)$ tidc (x) = ~(hic u + I 22 c 22 + /33C33), (29) 

with Iij = K n maf8ij /5. 

The dynamic equations of the system are found using the Euler-Lagrange equation: 

d fdL\ _ dL 
dr \dqi J %' 

where qi represents one of the dynamic variables {a±, a 2 , CI3, (ft, A}. We find 

2ir 5fciP 

ai = ai(Q 2 + A 2 ) - 2a 2 QA ai^ipH — cnai, (31) 

g n nK n p c ai 

a 2 = a 2 (0 2 + A 2 ) - 2ai^A a 2 ^ 2 pH L - £ - - c 22 a 2 , (32) 

g n nK n p c a 2 

2vr . 5/ciP c , , 

a 3 = a3^3pH 03303, (33) 

q n nK n p c a 3 



3 M r 2 + K 

J s = —K n m-^—^{al-a 2 2 )s\n[2{^-4>)}, (34) 
C = 0, (35) 



where we have defined 

2 

J, = /SI n n maia 2 A, (36) 

5 

2 

C = I A — -K n maia 2 Q, (37) 


and p = 3m/(47ri? 3 ), g ra = K n (l — n/5). The quantity J s is the spin angular momentum of the star 
and C is interpreted as the circulation in the locally nonrotating inertial frame. In the absence of 
viscosity, Kelvin's theorem ensures that the circulation is conserved. Note that the corotation case 
corresponds to A = 0. 

For a binary system in quasi-equilibrium, we require 

[W 

Oi = 0, <f> = V, and (ft = Q = * = J -j, (38) 

where the second equality comes from J s = and eq. (34), and the third equality follows from 
(11). It is of interest to note that in general relativity, quasi-equilibrium does not require the 
angular frequency of the figure, f2 = (ft, to be equal to the orbital angular frequency 0; rather, 
= (M/r 3 ) 1 / 2 , independent of the spin of the BH. 

The pressure terms in equations (31)-(33) can be conveniently rewritten in terms of the dynamic 
variables (LRS2): 



5ki P c m ( Ro\ 3 ^ n 
nn n p c q n Ro 



with R = ( aia2fl3 ) 1/3 



(39) 
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For n — > 0, we have k\jn — > 2/5 and K n — ► 1, so that 

5fei P c ^ 2^ 

Pc 



(40) 



For an incompressible ellipsoid, volume is conserved, R 3 = 010203 = R^; Differentiating this with 
respect to r gives ]Ci^i/ a j = 0. A second derivative yields X)i^i/ a i = Yl% ( a i/ a i) 2 ) which maybe 
combined with eqs. (31)-(33) to give 



2P r 



-1 



- E 



-2(tt 2 + A 2 )+2JM f^ + ^j +4vrp + ^ 



where we have used the relations + A2 + A3 = 2 and fen + fc 2 2 + &33 
quasi-equilibrium, we have dj = 0. 



(41) 

0. Of course, for 



3. Dimensionless Equations and Scaling Relations 

Before presenting our results on specific binary models, it is useful to cast the equations 
into a dimensionless form. This will shed light on how our results depend on various parameters 
(to, M, a, Rq) of the system. We define the following set of dimensionless variables: 

* = (42) 

a = "(i)" 2 - (43> 
A = A (i)" /2 - (44) 

(Tfi \ 
M) • (45) 

Since the tidal radius is expected to be of order Rq(M /to) 1 / 3 , we define the dimensionless orbital 
radius as 

r / to \ V 3 



fs wAM> ■ (46> 

In terms of these dimensionless variables, equations (31)-(33) can be written in the following form 
(for quasi-equilibrium) : 



3 



R y 1 (R \ 3/n 



= a 1 (n z + A z )-2a 2 nA--A 1 a 1 + ^ -cnfli, (47) 

l \ K J a\q n \ K J 

= a 2 (^ 2 + A 2 )-2d 1 M-^ 2 o 2 f§y + ^f§y /n -c 22 o 2 , (48) 

I \ K J a,2q n \ rl J 

R \ 3 1 (R \ 3/n 



= -Asa 3 (-) + - fn [ ir ) -3303, (49) 
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where R/Rq = (ai^a^) 1 / 3 , and 

h = f- 3 / 2 . (50) 

By making (40) and (41) dimensionless, the equations for the n = case may also be written in a 
dimensionless form. The dimensionless circulation is 



C 



(m 3 i? ) 1/2 5 



{a\ + a|)A - 2aia 2 £l . (51) 



In Newtonian theory, depends only on f, thus for given f and C, eqs. (47)-(51) completely 
determine the solution. However, this is not the case in general relativity since Cy depends not just 
on f. For example, 

*n = (" 2 " 3 ^ ) - ( 52 ) 

(The Newtonian case corresponds to if = 0). It is easy to see that relativistic effect enters through 
the combination K/r 2 , which can be expressed as 



K 



{r/M) 1/2 -a} 2 



r 2 (r/M) 2 -3(r/M) + 2a(r/M) 1/2 ' 



(53) 



where a = a/M. Therefore, in addition to the dimensionless Kerr parameter a, one could introduce 
another "parameter", r/M, which specifies how relativistic the orbit is. But of course this is not 
really a parameter, since r varies as the binary evolves. A more useful parameter is 



Ro = 



RojM/m) 1 / 3 _ Rq /mx2/3 

M m \m) ' 1 ' 



Roughly speaking, is the ratio of the tidal radius and the horizon radius of the BH. This 
parameter is convenient because it depends only on the mass and radius of the star and the mass 
of the BH. When Rq 3> 1, we expect that the binary will exist only in the Newtonian regime 
(i.e., the star will be disrupted before it enters into a relativistic orbit). But when Rq is small (of 
order unity), then the tidal disruption will be affected by general relativity. Note that r/M can be 
expressed in terms of f and ^ as 

£ = Rof. (55) 
Thus for given n, f, C, a and Rq, an equilibrium binary model is uniquely determined. 



4. Equilibrium Binary Models and Tidal Limits 



As discussed in §2, in the absence of viscosity, fluid circulation is conserved; gravitational 
radiation reaction, being a potential force, also preserves the circulation. In a coalescing binary, 
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we can determine the the relative importance of viscosity by comparing the orbital decay time 
td = \r/f \ and the tidal synchronization time t syD (Bildsten & Cutler 1992). The ratio is 

tJ^lf'ffflV, (56) 
t d \sm2aj \R J \m J v ; 

where a is the tidal lag angle, and is proportional to the viscosity. Kochanek (1992) and Bildsten 
& Cutler (1992) showed that t syn /td is greater than unity for coalescing NS-BH binary. Using 
microscopic viscosity for white dwarf matter (Nandkumar & Pethick 1984 and references therein) 
we can also show that the inequality t syn /td S> 1 is well satisfied for WD-BH binaries. Thus for 
the types of systems (NS-BH and WD-BH) of interest in this paper, the binary will likely evolve 
along a quasi-equilibrium sequence with constant C. We shall therefore focus on such models in 
this section. 



4.1. Roche- Riemann Model: Irrotational Equilibrium Sequences 

For the irrotational sequence, C = 0, the star has zero spin at large orbital radii (Q s = 0), and 
equation (51) implies A = 2a\a2^l/{a\ + a|). In Figure 2, four equilibrium sequences are depicted 
illustrating the effect of general relativity. Note that it is most convenient to parameterize the 
sequence in terms of 02/01 since r is double- valued. We see that the star becomes more distorted 
as f decreases. At the same time, its central density p c decreases. This implies that the tidal 
field stablizes the star against collapse. Physically, this stablization arises from that fact that the 
self-gravity of a distorted object is weaker, and the star expands in order to satisfy hydrostatic 
equilibrium (see, e.g., Lai 1996; LRS1; Brady & Hughes 1997; Wiseman 1997; Baumgarte et 
al. 1998; Thorne 1998; Flanagan 1998; Bonazzola et al. 1999a,b). 

For each sequence, there exists a turning point, corresponding to the tidal limit, r t id e , below 
which no equilibrium solution is possible. From Fig. 2, we see that stars with finite Ro's are tidally 
disrupted at a larger f than purely Newtonian stars (with Rq — ► 00), an incarnation of the old adage 
that GR means "stronger" gravity. We also see that the tidal limit fy^e is smaller for positive finite 
d (i.e., the orbit is prograde with respect to the BH spin), while f t id e is larger for more negative d. 
This effect of the Kerr parameter d might also be guessed from the effect it has on the innermost 
stable circular orbit (ISCO): The ISCO is largest for a = — 1, 

^isco — 9i\^f, and smallest for d — 1, 
risco = M. The same qualitative effect applies in tidal disruption; the strongest tidal gravity is seen 
for the a = — 1 case. 

The dependence of the tidal limit on d can be understood by examining the tidal field (eqs. [18]- 
[21]). The spin parameter d enters the tidal field only through K/r 2 . Figure 3 shows contours of 
constant K/r 2 in the d — (r/M) plane. Note that only the region outside the ISCO (r > rj sco ) 
is of physical interest. At the ISCO, it can be shown that K/r 2 = 1/3, independent of d (This 
coincidence was first noted by Fishbone 1973). We see from Fig. 3 that at a given r, the quantity 
K/r 2 is smaller for larger d, implying that the tidal field is weaker for larger, positive d. 
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4.2. Irrotational Tidal Limit 

As discussed above, the turning point in the equilibrium sequence is of special significance: it 
corresponds to the tidal limit (also called "Roche-Riemann limit"), ftide- To determine this limit, 
we require 8f / 8(0,2/ '01) = 0, which adds another equation to be solved together with eqs. (47)-(49). 

We have calculated f t id c as a function of Rq for a series of polytropic indices, n = {0, 0.5, 1, 1.5} 
and three Kerr parameters, a = { — 1, 0, 1}. The results are are shown in Fig. 4 and Fig. 5. Selected 
numerical values are also listed in Table 1. For given n and a, there is a critical Rq = -Ro.crit, for 
which Hsco = r t ide> i-e., the inner-most stable circular orbit (ISCO) coincides with the tidal radius. 
For Rq < i?o,crit> the binary encounters ISCO before tidal disruption; In this case, tidal disruption 
can still occur outside the BH's event horizon, but it takes place on a plunging orbit rather than on 
a circular orbit as in the Rq > -Ro,crit case. Our results in Table 1 and Figs. 4-5 therefore have been 
truncated at Rq = ^o,crit- Note that for R = ^o,crit( a ) ( so that tidal disruption occurs at ISCO), 
the tidal limit f t id e is independent of a. This is because in the framework of our calculations, general 
relativistic effects come into the equations only through K/r 2 , which is equal to 1/3 (independent 
of a) when evaluated at r; sco (a) (see Fig. 3). 

Figure 6 depicts the central density p c and the distortion (a = 02/01 or 03/01) of the star at 
the tidal limit as a function of Rq. Only the n = 1 results are shown as an example. We see that 
Pc,tide and a t ide (evaluated at r = rtide) are rather insensitive to Ro and d, although the values of 
ftidc for different Rq and a can differ by as much as 10% (see Table 1 and Figs. 4-5). 



4.3. The Effect of Initial Spin 

When the star has an intrinsic spin f2 s at large orbital radii, its (conserved) circulation is finite. 
Since at r — > 00, the star is axisymmetric (ai = 02), and the figure rotation Q — > 0, we have (see 
eq. [51]) 

C = --KnO^s, (57) 

where we have identified Q s = — A(r — > 00). Note that when fl s is positive (i.e., the stellar spin is 
in the same direction as the orbital angular momentum), C is negative. The maximum spin that 
a uniformly rotating star can sustain without shedding mass from its equator is given by (e.g., 
Friedman, Ipser & Parker 1986; Cook, Shapiro Sz Teukolsky 1992) 

(m/i?o) ' 

For a given Cl s , the dimensionless circulation C can be calculated by using the equations for com- 
pressible Maclaurin spheroid (LRS1). 

To illustrate the dependence of the tidal limit on the initial stellar spin, Figure 7 shows ftide 
as a function ^0 for n = 0.5, d = and several different values of Cl s . For negative Cl s (positive C), 
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the internal rotation A is larger (compared to the C = case), and this induces larger distortion 
to the star, and thus the star gets disrupted at a larger radius. For small, positive Cl s , the internal 
rotation A is reduced, and thus ftide is smaller. But as f2 s becomes larger (Cl s <; 0.3), A also gets 
larger and ftide is increased over the Cl s = value. 



4.4. Roche Model: Corotating Binaries 

Although it is likely that equilibrium sequences with constant circulation should describe the 
quasi-static evolutionary phase of NS-BH binaries (see the beginning of §4 for discussion), some 
degree of spin-orbit synchronization can still be achieved at small orbital radii through an anoma- 
lously high (turbulent) viscosity (so that sin 2a is not much less than unity). This is especially true 
for WD-BH binaries, for which corotating models (with A = 0) may be relevant. For completeness, 
Table 2 lists selected values of the tidal limit (of "Roche limit") as a function of Rq and a for 
polytropic index n = and n = 1.5. For the incompressible model (n = 0), our result agrees with 
that of Fishbone (1973) (although a different parametrization was adopted in Fishbone's paper). 
Similar to the irrotational models, we find that the GR effect makes ftide larger than the Newtonian 
value, and that ftide is larger for a = — 1 than for a = 1. 



5. Tidal Effect on the Orbit 

The previous sections are based on the assumption that the star m moves in a geodesic of the 
BH metric. In reality, there are two corrections to the geodesic. First, even if m were a point mass, 
its finite mass can modify the spacetime of the system. For r>M|m, the problem can be solved 
by post-Newtonian expansion (e.g., Blanchet et al. 1998 and references therein) or BH perturbation 
theory (e.g., Mino et al. 1997); when r is comparable to M + m, nonperturbative techniques are 
needed (see Kidder, Will & Wiseman 1993, Buonanno & Damour 1998 and references therein) and 
definitive result will only come from numerical relativity. Second, there is a correction due to the 
finite size of the star. Here we focus on the second correction, i.e., the tidal effect on the orbit. 



The tidal interaction energy is given by 



2r 3 

(see eq. [29]). The dimensionless ratio 

r Wtide K« /m\2/3 1 



K 



(59) 



(60) 



measures the correction of the tidal effect on the orbital motion (deviation from geodesic) 2 . Since 



2 In Newtonian theory, the orbital angular freqency is given by fi = \{M + m)/r 3 ] 1 ^ 2 (l + 35/2). 
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(2a 2 — a\ — a 2 ) and (a 2 — S§) are of order 1/f 3 , and K/r 2 is of order M/r < 1, we see that 5 scales 
as 1/f 5 . Note that even near the tidal radius (f ~ 1), 5 can be quite small if m/M is small. 

In the absence of tidal coupling, the dynamical instability limit r& yn of the orbit is located at 
ISCO. We assume that ri SC0 is given by the standard M m result. (When m/M is finite, see 
Buonanno &: Damour 1999 and references therein.) Tidal coupling is expected to shift rdyn outward 
compared to nsco (see Lai <fc Wiseman 1996). To estimate the correction, (r& yn — r- lsco ), we write 
the total equilibrium energy of the binary system as 

S = mE + AS, (61) 

where E is the energy per unit mass of a test particle orbiting a Kerr BH (see eq. [9]), and AS is 
the correction due to the finite size of the star. Minimizing S would give r^ yn . Unfortunately there 
is no unambiguous way of writing down the total energy S in general relativity. However, we can 
still attempt to determine how (rd yn — Hsco) scales with various parameters of the system. From 
Newtonian consideration, AS is certainly of order (— Wtide) — The negative sign arises because the 
tidal attraction increases the kinetic energy of the orbital motion, which more than compensates 
for the neative tidal potential energy. We therefore set AS = — (3W t ide, with f3 of order unity. 
Assuming that r^ yn is close to ri sco , we find 3 : 

= r dyn - Hsco _ £ 2 ( ™ \ 2 / 3 



M 



C(iK n Rl{^) (2a\-a 2 2 -al). sco , (62) 



where (2S 2 — a?, — S3) is evaluated at r = ri sco . The constant C depends on r\ sco /M and d 2 E/dr 2 
(evaluated at r isco ): C = 0.116 for a = 0.998, and 0.106 for a = and 0.077 for S = -1. Although 
we cannot determine the precise value of (3, we expect the scaling in eq. (62) to be valid in general. 

The biggest correction to r^ yn occurs when tidal disruption occurs just inside the ISCO. For 
example, when Rq = 2, a = and n = 0, we have e = 0.53/3(ra/M) 2 / 3 (assuming m/M is small). 



6. Possible Applications to Gamma-Ray Bursts 

In this section we discuss possible applications of our results to Gamma- Ray Bursts (GRBs). 



6.1. Black Hole — Neutron Star Binaries 

Coalescence of BH-NS binaries has frequently been invoked as the central engine of GRBs 
(e.g., Narayan, Paczynski & Piran 1992; Meszaros, Rees & Wijers 1998; Piran 1998). It is expected 



3 In taking derivative of Wtidc with respect to r, we have assumed that K/r 2 is roughly constant, which we have 
checked to be reasonable as long as a is not too close to 1. Also we have found that (2a? — a? 2 — a§) is approximately 
equal to twice of (of — a§), and both scale approximately as r~ 3 . Relaxing these approximations would give rise to 
a factor of order unity which we can absorb into the definition of f3. 
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that a thick disk or torus will result from the disruption of the NS. To account for the temporal 
properties of observed GRBs, the disk must be long-lived (compared to the orbital timescale) (e.g., 
Meszaros, Rees & Wijers 1998; van Putten 1999), and therefore the tidal disruption must occur 
outside the ISCO, i.e., r t id e > nsco- Approximating the NS as a n = 0.5 polytrope, we find that 
the tidal radius for critical R = i?o,crit (where r ti d c = r- lsco ) is f t id e = 2.56 (see Table 1). Therefore 
the condition r t id c > nsco translates to 

M (M\ /i? \ 3/2 /6M \ 3/2 , , 

For canonical neutron star parameters (m = 1.4M , Rq = 10 km, corresponding to Ro/m ~ 4.84), 
tidal disruption occurs outside ISCO only if M ^ 4.4M (for a = BH) or M <> 64M (for a = 1 
BH). 

We can calculate the maximum energy available in the disk following a disruption by comparing 
the binding energies at r ti( jc and at rj sco . Assuming that the disk material follows geodesic, the 
maximum energy release per unit mass of accreting gas is given by 

AE = E(r tide ) - E(r isco ), (64) 

where E(r) is given by eq. (9). Note that AE depends only on Rq. Figure 8 depicts AE as a 
function of the "effective" mass ratio (M/m)(R /5m)- 3 / 2 = (5/R ) 3 / 2 for different value of BH 
spin parameter a (a = —1, 0, 0.5, 0.9, 0.998 and a — ► 1; Note that a = 0.998 corresponds to the 
maximum value of a that can be achieved in accretion; Thorne 1974). Clearly, AE = when 
M/m = (M/m) cr it since r t id c = Hsco- When (M/m) <C (M/m) cr i t , AE approaches the standard 
radiation efficiency for thin disk, AE = l — E(ri sco ). The dependence of AE on a mainly comes from 
Hsco since rtide is relatively insensitive to a (see Figs. 4-5 and Table 1). Note that our calculation is 
valid only when M/m 3> 1, so the results for M/m less than a few should be considered suggestive 
only. 

Recent observations indicate that some GRBs are highly energetic, with the apparent emitted 
energy in gamma rays reaching 10 53 erg (e.g., Kulkarni et al. 1999 and references therein). Un- 
certainty in beaming precludes accurate determination of the true gamma ray energy. It has been 
suggested that the total energy in the fireball may be as high as 10 54 erg (e.g., Kumar 1999). If 
this number holds up, we may rule out NS-BH binary coalescence as the central engine of GRBs 
(see Fig. 8). 



6.2. Black Hole - White Dwarf Binaries 

For a BH-WD binary, tidal disruption occurs outside ISCO (r ti d c > Hsco) when 

M (M\ ,( R \ 3/2 /6M\ 3/2 , x 

™ n™)*,- 2 ' 18xl0 U4) UJ ■ < 65 > 
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where we have adopted the n = 1.5 polytropic model for the WD, and used f t id c = 2.173 (see Table 
1) for R = Rq. cr it. Although there is no evidence for massive BHs with mass around 1O 5 M (e.g., 
Richstone et al. 1998), the observational constraint on such BHs is also weak. If they do exist, they 
may capture WDs and the merger of the BH-WD binary may lead to a GRB (R. Blandford 1997, 
private communication). For (M/m) <C (M/m) crit , a Newtonian description for tidal disruption is 
adequate; but when M/m is close to (M/m) cr i t , a general relativistic treatment is necessary. Since 
M>mis well satisfied, the model presented in this paper should provide an accurate description 
of BH-WD binaries. 

Figure 8 shows the maximum energy available in a BH-WD tidal disruption. Unlike BH-NS 
binaries, the limiting value of AE, corresponding to the thin disk limit, AE = 1 — E(ri sco ), can be 
easily achieved for a large range of mass ratios. 

7. Discussion 

Quasi-equilibrium models of neutron star binaries provide an accurate description of the late 
stage of inspiralling binary stars prior to the dynamical merger. While quasi-equilibrium NS-NS 
binary models (the so-called "Darwin- Riemann problems") have been successfully constructed in 
general relativity (see Bonazzola et al. 1999a,b for a review), only Newtonian models of NS-BH 
binaries have been studied so far (see Uryu & Eriguchi 1999). The semi-analytic Roche- Riemann 
models studied in this paper provide useful insight on the relativistic tidal effects in coalescing 
NS-BH and WD-BH binaries. 

We note that our model is exact when the following conditions are satisfied: (i) M S> m [so 
that the spacetime is given by the Kerr metric of the BH (§2), and the tidal correction to the orbit 
is negligible (§5)]; (ii) Ro/m 3> 1 (so that we can treat the fluid star as a Newtonian object); (iii) 
The radius of curvature (r 3 /M) 1 / 2 » Rq [so that high-order tidal effect can be neglected (§2); 
At the r = r ti< j c ~ (M/m) 1 / 3 ^, this condition translates to (Ro/m) 1 / 2 3> 1, and is equivalent to 
condition (ii)]; (iv) n = [The ellipsoid ansatz is exact when the external force is a linear function of 
spatial coordinates, as required by (iii), and when the fluid is incompressible (Chandrasekhar 1969)]. 
Clearly, for WD-BH binaries, conditions (i)-(iii) are well-satisfied, and the ellipsoid ansatz is a good 
approximation for polytrope of n = 1.5 (appropriate for low-mass WD) (LRS1 4 ). Therefore our 
findings that general relativitic effect increases r t id e and that r t jd e is largest for a = — 1 (regtrograde 
BH) are expected to hold. 

For NS-BH binaries, conditions (i)-(iii) are satisfied only approximately, we must apply our 
results with caution. In particular, post-Newtonian correction to the self-gravity is likely to be im- 
portant. The height H of tidal bulge on the star can be estimated from (MRo/r 3 ) [l + 0(K/r 2 )] ~ 



4 Note that for very compressible fluid (large n), the tidal limit corresponds to the point where mass transfer starts; 
see Uryu & Eriguchi (1999). 
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(mH/Rl) [1 + O(m/R )]. Setting H ~ i? gives the tidal radius 

« ^\ 
rtide ~ ito — 
\m J 

While relativistic tidal field (the if/r 2 term) increases rtide ; the post-Newtonian self-gravity effect 
(the m/Ro term) tends to reduce rtide- Since m/Ro for a NS is of order 20%, comparable to 
K/r 2 (which is always less than 1/3; see Fig. 3), the post-Newtonian self-gravity effect cannot be 
neglected. Note that our finding that r tic j c is largest for a = —1 is still expected to be valid. 

P.W. acknowledges support from the NSF REU program at Cornell, and D.L. thanks the 
Alfred P. Sloan Foundation for a research fellowship. 
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Table 1. Tidal Limits for Irrotational Binaries 
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Table 1 — Continued 
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Note. R - %m 2/3 

(m, i?o are the mass and radius of 
the star, M is the BH mass), r t id e = 
(m) 1 ^' " is tne Kerr parame- 
ter of the BH, and n is the poly- 
tropic index. The last item in each 
column corresponds to Rq = i?o,crit 
at which r ti d c = Hsco or Rohide = 
(jwOisco - 
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Table 2. Tidal Limits for Corotating Binaries 
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Fig. 1. — Sketch of the Roche- Riemann binary model. The star is represented by an ellipsoid with 
mass m. The BH has mass M and Kerr parameter a (The BH spin is aligned with the orbital 
angular momentum). The angle * is shown with quotes since the drawing is only correct in the 
Newtonian limit (see text). 



-23- 




Fig. 2. — Equilibrium sequences for irrotational binaries (C = 0). The axis ratio a<ija\ and the 
central density p c [in units of its value at infinite binary separation, po = Pc(r — > oo)] are plotted 
against the dimensionless orbital separation f = (r / Ro)(m/M) 1 l' i . The polytropic index is n = 1. 
The dashed curves represent the Newtonian result, with Rq — > oo. The three solid curves include 
the effects of general relativity (with a = 1, 0, — 1 from left to right), all for Rq = 4. The turning 
point of each curve, with f reaching a minimum, is marked by an asterisk, and corresponds to the 
tidal limit. Note that the values of 02/01 and p c /po at f = ftide are similar for different sequences, 
but they are not identical (see Fig. 6). 
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a 



Fig. 3. — Contours of constant K/r 2 in the a — (r/M) plane. This quantity determines the 
magnitude of relativistic effect on the tidal distortion. The number on each curve gives the value of 
K/r 2 . For physically reasonable r/M, the quantity K/r 2 is positive and decreases with increasing 
a and r/M. Note that at the innermost stable circular orbit (ISCO), K/r 2 = 1/3 is independent 
of a. Thus the heavy solid line also gives the value of r; SC o as a function of a. The pole of K/r 2 
(where it becomes infinite) is hidden inside the ISCO except when a = 1. 
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2 ^ i i i i—i i ' i i 

'0 1 2 3 4 5 6 7 

Fig. 4.— Tidal limit f tidc = (r tidc / R )(m/M) 1 / 3 as a function of R = (R /m)(m/M) 2 / 3 for 
irrotational binaries (C = 0). The two sets of solid curves correspond to polytropic index n = 
and n = 1 respectively. Results for three values of Kerr parameter (a = —1,0, 1) are shown. Each 
solid curve is truncated at Rq — ^o,crit, where — ?*isco [The short-dashed lines represent the 
ISCO for different Kerr parameters, i.e., R^r = (r/M); sco ]. The two horizontal dotted lines show 
the limiting rude a t the truncation. Two additional long-dashed lines show the Newtonian limit of 
the tidal limit. For n = 0, the Newtonian limit is rude = 2.502 and the ISCO limit is rtide = 2.757. 
For n = l, the Newtonian limit is f t id e = 2.140 and the ISCO limit is r t ide = 2.366. 
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Fig. 5. — Same as Fig. 4 except for n = 
ftide = 2.320 and the IS CO limit is ftide = 
and the IS CO limit is ftide = 2.173. 



0.5 and n = 1.5. For n = 0.5, the Newtonian limit is 
2.560. For n = 1.5, the Newtonian limit is f t jd e = 1.960 
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Fig. 6. — The central density [in units of po = p c (r — > oo)] and the axis ratio (a = 02/01 or 03/01) 
of the star at the tidal limit (f = "r t ^ e ) as a function of Rq for a = {1,0,— 1}, n = 1 and C = 0. The 
solid curves in the lower plot are for a = 02/01 and the dashed curves are for a = 03/01. From left 
to right, the curves are a = {1,0,-1}. The horizontal dotted lines correspond to the Newtonian 
limit. 
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Fig. 7. — Tidal limits in models with initial rotation (fl s ^ 0) for n = 0.5 and a = 0. As in Figs. 4-5, 
each curve is truncated at Rq = -Ro,crit> where r t id e = Hsco [The short-dashed line represents the 
ISCO, i.e., fRo = (r/M) isco }. 
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Fig. 8. — The maximum energy available (per unit mass of accreting matter), AE = E(r t ^ e ) — 
-E(nsco)) in the accretion disk following the tidal disruption, is plotted against the "effective" mass 
ratio. The numerical value adjacent to each solid curve indicates the Kerr parameter a. (a = 
0.998 is the Thorne limit for Kerr BHs fed by accretion disks.) The dotted horizontal lines indicate 
the limiting energy for disruption at oo for each of the Kerr parameters. For BH-NS binaries, the 
bottom horizontal axis shows the mass ratio scaled for Rq/tti = 5 (Results for M/m less than a few 
can only be viewed suggestively, and results below M/m = 1 is not physically meaningful). For 
BH-WD binaries, the top horizontal axis shows the mass ratio scaled for Rq/tti = 10 4 . The label on 
the right gives the total energy in ergs for m = 1.4M©. Note that in this figure, r tic j is calculated 
for n = 0, but the dependence of AE on n is weak. 



